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Using a framework in which all elements are constrained by Dyson-Schwinger equation studies 
in QCD, and therefore incorporates a consistent, direct and simultaneous description of light- and 
heavy-quarks and the states they constitute, we analyze the accuracy of 5'?7(4)-flavor symmetry re- 
lations between Trpvr, KpK and DpD couplings. Such relations are widely used in phenomenological 
analyses of the interactions between matter and charmed mesons. We find that whilst 5'f/(3)-flavor 
symmetry is accurate to 20%, 5(7(4) relations underestimate the DpD coupling by a factor of five. 

PACS numbers: 14.40.Lb, ll.15.Tk, 12.39.Ki, 24.85.-|-p 



I. Introduction. Hadrons in-medium arc the focus of 
intense theoretical and experimental activity. The chief 
motivation in heavy-ion collisions is a better understand- 
ing of QCD's deconfined phase, viz. the putative quark- 
gluon plasma, its chiral restoration phase transition and 
associated order parameters. Whilst an enhancement of 
charm and strangeness in the quark-gluon phase is pre- 
dicted to lead to the copious production of mesons 
[H at the large hadron collider, J/ip suppression has long 
been suggested as an unambiguous signature for quark- 
gluon plasma formation Notwithstanding ongoing 
debates about charmonia production mechanisms and a 
wide range of suppression effects, much effort is sensi- 
bly dedicated to understanding the complicated final- 
state interactions which occur after hadronization of the 
plasma; see, e.g., Ref. Q. 

Charmed-meson interactions with nuclear matter will 
also be studied at the future Facility for Antiproton and 
Ion Research (FAIR) and possibly at Jefferson Labora- 
tory (JLab). Low-momentum charmonia, such as J/ip 
and and D'*-' mesons can be produced by annihilation 
of antiprotons on nuclei (FAIR) or by scattering electrons 
from nuclei (JLab). Since charmonia do not share valence 
quarks in common with the surrounding nuclear medium, 
proposed interaction mechanisms include: QCD van der 
Waals forces, arising from the exchange of two or more 
gluons between color-singlet states 3; and intermediate 
charmed hadron states |5|, such that hadronic 
vacuum polarization components of the J/ip interact with 
the medium via meson exchanges . 

A kindred approach is applied to low-energy inter- 
actions of open-charm mesons with nuclei, which may 
create a path to the production of charmed nuclear 
bound states (D-mesic nuclei) These studies 

rely on model Lagrangians, within which effective inter- 
actions are expressed through couplings between D^*^- 
and light-pseudoscalar- and vector-mesons. The models 
are typically an SU{4:) extension of light-flavor chirally- 
symmetric Lagrangians. Most recently, exotic states 



formed by heavy mesons and a nucleon were investigated, 
based upon heavy- meson chiral perturbation theory [ld| ■ 
In that study a universal coupling, g^^, between a heavy 
quark and a light pseudoscalar or vector meson was in- 
ferred from the strong decay D* — )- Dn, cf. Ref. [lT| . 

In the context of chiral Lagrangians, it is natural to 
question the reliability of couplings based on SU{A) sym- 
metry. Flavor breaking effects are already known to oc- 
cur in the strange sector and should only be expected to 
increase when including charm quarks. The order of mag- 
nitude of this larger symmetry breaking is signalled by 
the compilation of charmed couplings in Ref. |3| , where 
SU{A) relations are shown to be violated at various de- 
grees (ranging from 7% to 70%) in couplings between two 
heavy mesons and one light meson. No states containing 
a s-quark were considered. 

Herein, we study a different quantitative measure, 
based upon ratios between the DpD, KpK and Trpir cou- 
plings; namely, a difference between the same coupling in- 
volving either a c-, s- or light-quark. We are motivated by 
the notion that the KpK and DpD systems are dynami- 
cally equivalent in the sense that the heavier quark acts as 
a spectator and contributes predominantly to the static 
properties of the mesons, whereas the exchange dynam- 
ics is mediated by the light quarks. In practice, the sym- 
metry idea is expressed by implementing gopD — 9 KpK 
in the meson-exchange models d, Q. The Trpn coupling 
provides a well-constrained benchmark. 
II. DSE Pramev^rork. Our primary object of interest 
is a phenomenological coupling that relates the transition 
amplitude of an initial pseudoscalar H — Qf -meson, Q — 
c, s and f = u,d, to an identical meson via emission of 
an off-shell p. The matrix element for this transition is 

{H{p2)\p{P,X)\H{p,)) =gHpH ex-P, (1) 

an expression which defines the dimensionless coupling of 
the two pseudoscalar mesons to a vector meson with mo- 
mentum P = P2 —pi and polarization state A. The decay 
p — > TTTT is also described by such a matrix element. How- 
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ever, there is no associated physical process when < 



Amjj and Pi = P2 = (N.B. A Euchdean metric 

is used: {7^,7,.} = 2(5^^; 7^=7^; a-6 = SLi'**^*; 
and tr[757^7i.7p7(T] = -4e^^po-, £1234 = 1- For a space- 
Uke vector P^ , > 0.) Nevertheless, a coupling of this 
sort is employed in defining p-meson-mediated exchange- 
interactions between a nucleon and pseudoscalar strange- 
or charm-mesons. In such applications: the off-shell p- 
meson's momentum is necessarily spacelike; and a cou- 
pling and form factor may be defined once one settles on 
a definition of the off-shell p-meson. 

Symmetry-preserving models built upon predictions 
of QCD's Dyson-Schwinger equations (DSEs) provide a 
sound framework within which to examine heavy-meson 
observables pH - fTsj . Such studies describe quark propa- 
gation via fully dressed Schwinger functions, which has a 
material impact on light-quark characteristics [l6j . 

At leading-order in a systematic, symmetry-preserving 
truncation scheme [ItI , one may express Eq. ([T|) as 



9HpH e ■ P = trcD 



d*k 



Tnik; ki)SQ{kQ) 



(2^)4 

X fHik;-k2)Sfik'f)e^*-fp{k;-P)Sfikf) , (2) 

where S represent dressed-quark propagators for the in- 
dicated flavor and are meson Bethe-Salpeter ampli- 
tudes (BSAs), with H = tt,K,D. In Eq. ©: the trace 
is over color and spinor indices; kq — k + wipi^k'j = 
k + wipi — P2, kf = k — 'W2P1, where the relative- mo- 
mentum partitioning parameters satisfy 1^1+102 = 1; and 

is the vector-meson polarization four- vector. This ap- 
proximation has been employed successfully; see, for in- 
stance, applications in Refs. |T^. HgI. [l8l - l22l | . 

We simultaneously calculate the D-, K- and p- meson 
leptonic decay constants via [l3j : 



Pp/ff = trcD 
1 



757MXH(fc;P) , 



(3) 



Mpfp 



:trcD 



(2^)4 

d k „ , , , , 

^7.x^(fc;i^), (4) 



where x(fc; P) = Sf^ {k + wiP)T{k; P)Sf^ {k - waP). The 
BSAs are canonically normalized; viz.^ for pseudoscalars 



2P^ = 



d 



OK, 



■n(p,A'; 



K=P 



(5) 



n(P,if) =trcD 



d'^k 



fo-{k;^P)Sf,ik + wiK) 



(27r)4 

xro-{k;P)SfAk-W2K), (6) 

with an analogous expression for the p [l^. 

The solution of QCD's gap equation is the dressed- 
quark propagator, which has the general form 

S{p) = ^ij.pav{p'^) + <Ts{p^) = l/[i-fpA{p^) + B{p^)] . 

(7) 

For light-quarks, it is a longstanding DSE predic- 
tion that both the wave-function renormalization. 



^{p^) ~ and dressed-quark mass-function, 

Af(p^) = B{p'^)/A{p'^) = crs{p'^)/avip'^), receive strong 
momentum-dependent modifications at infrared mo- 
menta: Z{p^) is suppressed and AI{p'^) enhanced. These 
features are characteristic of dynamical chiral symmetry 
breaking (DCSB) and, plausibly, of confinement. (N.B. 
Eqs. ([5]), ^ represent the quark propagator S{p) as an 
entire function, which entails the absence of a Lehmann 
representation and is a sufficient condition for confine- 
ment [13, [l^l-) The significance of this infrared dressing 
has long been emphasized [Tsj : e.g., it is intimately con- 
nected with the appearance of Goldstone modes [la| . The 
predicted behavior of Z{p^), M^p^) has been confirmed in 
numerical simulations of lattice-regularized QCD [l^jHE]- 
Whilst numerical solutions of the quark DSE are read- 
ily obtained, the utility of an algebraic form for S'(p), 
when calculations require the evaluation of numerous in- 
tegrals, is self-evident. An efficacious parametrization, 
exhibi ting the aforementioned features and used exten- 
sively dl, [13, [H, [2^ , is expressed via 



as{x) 2 m J"(2(a; + ffi^)) 

+T{bix) T{b3x) [bo + b2J'{ex)] 
1 



(Tvix) 



[1 - Ti2{x + rh'^))] 



(8) 
(9) 



m/A, J-{x) 

\2 . 



^<^s{p^) and ayix) = av{p^)- The parame- 



[1 
o2^ 



exp(-a;)]/x, 



with X 

a six) = Aasyp-) ana ay I 

ter values were fixed [l3| by requiring a least-squares fit 
to a wide range of light- and heavy-meson observables, 
and take the values: 



m/ 



■^0 



"2 



b^ 
"3 



u^d 0.00948 0.131 2.94 0.733 0.185 
s 0.210 0.105 3.18 0.858 0.185 



(10) 



At a scale A — 0.566 GeV, the current-quark masses 
take the values m„ = 5.4 MeV and ms — 119 MeV, and 
one obtains the following Euclidean constituent-quark 
masses [30|: = 0.36 GeV and Mf = 0.49 GeV. (N.B. 
e = 10^"* in Eq. dH) acts only to decouple the large- and 
intermediate-^^ domains |18l|.) 

We note that studies which do not or cannot imple- 
ment light-quark dressing in this QCD-consistent man- 
ner invariably encounter problems arising from the need 
to employ large constituent-quark masses and the associ- 
ated poles in the light-quark propagators [27| . This typ- 
ically translates into considerable model sensitivity for 
computed observables [l5| . 

Whereas the impact of DCSB on light-quark propa- 
gators is significant, the effect diminishes with increas- 
ing current-quark mass (see, e.g.. Fig. 1 in Ref. [isj). 
This can be explicated by considering the dimensionless 
and renormalization-group- invariant ratio <^f := (jf/Mj^ 
where cr/ is a constituent-quark cr-term: measures the 
effect of explicit chiral symmetry breaking on the dressed- 
quark mass-function compared with the sum of the effects 
of explicit and dynamical chiral symmetry breaking. Cal- 
culation reveals [2^: <;„ = 0.02, q = 0.23, = 0.65, 
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c^b = 0.8. Plainly, <;f vanishes in the chiral limit and re- 
mains small for light quarks, since the magnitude of their 
constituent mass owes primarily to DCSB. On the other 
hand, for heavy quarks, — >■ 1 because explicit chiral 
symmetry breaking is the dominant source of their mass. 
Notwithstanding this, confinement remains important for 
the heavy-quarks. These considerations are balanced in 
the following simple form for the c-quark propagator: 



Sc{k) = 



'ij ■ k + Mc 



He /Ml) . 



(11) 



which implements confinement but produces a 
momentum-independent c-quark mass-functiori; namely, 
cr^(fc2)/cr|(fc2) = Mc. We use Mc = 1.32 GeV 

A meson is described by the amplitude obtained from 
a homogeneous Bethe-Salpeter equation. In solving that 
equation the simultaneous solution of the gap equation 
is required. Since we have already chosen to simplify the 
calculations by parametrizing S\p), we follow Refs. [TTI . 
[Tsl - flll l and also employ that expedient with Th{p)- 

In this connection, the quark-level Goldberger- 
Trciman relations derived in Ref. psl ] motivate and sup- 
port the following paramctrization of the tt and K BSAs: 



(12) 



where Bt^ k '■— ^ulm ->o ^^'l obtained from 
Eqs. - ® through the replacements 5^ 5J = 0.204, 
^0 ~^ ^0 ^ 0.319, which yield computed values = 
146 MeV, fx = 178 MeV Equation ^ expresses 

the fact that the dominant invariant function in a pseu- 
doscalar meson's BSA is closely related to the scalar piece 
of the dressed-quark self energy owing to the axial- vector 
Ward-Takahashi identity and DCSB. 

Regarding the p meson, DSE studies [13, [11] indicate 
that, in applications such as ours, one may effectively use 



1_ \ P' , (13) 



r(^(fc;P)- 7 



p2 



namely, a function whose support is greatest in the in- 
frared. Similarly, for the D meson we choose: 



exp(-fcVt^i,) 



Ml 



D 



(14) 



The normalizations, A/'p, Md, are obtained from 
Eqs. (O, (IH) and simultaneous calculation of the weak 
decay constant in Eqs. ([3]), (HJ. In the expression for the 
coupling, Eq. ([T|), as well as in Eqs. ©-([l]), we follow the 
momentum-partitioning prescription of Ref. [llj, which 
leads to w1 = 0.79; viz., most but not all the heavy- 
light-meson's momentum is carried by the c-quark. We 
note that Poincare covariance is a hallmark of the direct 
application of DSEs to the calculation of hadron prop- 
erties. In such an approach, no physical observable can 



depend on the choice of momentum partitioning. How- 
ever, that feature is compromised if, as herein, one does 
not retain the complete structure of hadron bound-state 
amplitudes js^l ■ Any sensitivity to the partitioning is an 
artifact arising from our simplifications [111, [l^l ■ 
III. Results. The D-meson's width parameter is de- 
termined via analysis of relevant leptonic and strong 
decays: wr, = 1.63±0.10GeV for mo ^ 1.865GcV 
yields fo = 206 ± 9 MeV [H and go-Dn = IS.ltW cf. 
17.9 ± 1.9 HH- For the p, we use _Wp = 0.56 ± 0.01 GeV 
and 



0.38, both determined [ij] via a least-squares 
fit to an array of light-light- and heavy-light-meson ob- 
servables with rup = 0.77 GeV. Using Eqs. 0, (0 and 
([6|) , one therewith obtains fp ~ 209 MeV, cf. experiment 
216 MeV, which follows from the e+e" decay width j33l |. 

With the width parameters fixed, we computed the 
DpD, K pK and ttjOtt couplings in impulse approxima- 
tion, following Eq. Our results are depicted in Fig.[T] 
Notably, we compute the amplitude directly: at all values 
of and current-quark mass. We do not need to resort 
to extrapolations, neither from spacelike timelike mo- 
menta nor in current-quark mass, expedients which are 
necessary in some other approaches [1, . 

The behavior of q^„^(P-^) provides a context for our 



results. Experimentally 



r(— TOp) ~ 6.0; and 



: g-npiT\ ■■"pj 

the best numerically-intensive DSE computation avail- 
able produces [l^] 57rp7r(— "ip) = 5.2. Our algebraically- 
simplified framework produces .97rp7r(— w^) = 4.8, just 8% 
smaller than the latter, and a P^-dependence for the cou- 
pling which closely resembles that in Ref. [11]; e.g., both 
are smooth, monotonically decreasing functions and our 
value of (7^p7r(— TOp)/.g7rp7r("T-p) = 0.14 is just 10% smaller. 
On the domain P^ e [- 



1.84- 1.45s 
- 0.75s + 0.085s2 



(15) 



provides an accurate interpolation of our result. If one 
insists on a monopole paramctrization at spacelike- P^, 
then a monopole mass of A^p^ = 0.61 GeV provides a fit 
with relative-error-standard-deviation = 5%. 

In the case of exact SU{3) symmetry, one would have 
gKpK = 37rp7r/2. It is clcar from the figure that the 
assumption provides a fair approximation to our result 
on a domain which one can reasonably consider as rele- 
vant to meson-exchange model phenomenology; viz., on 



P2 G [- 



g] the error ranges from (—10)^40%. On 



this domain an accurate interpolation is provided by 



9Kpk{s) 



0.94 - 0.62s 
1 + 0.55s -0.16s2' 



(16) 



If one insists on a monopole paramctrization at spacelike- 
P^ , then a monopole mass of t^xpK = 0.77 GeV provides 
a fit with relative-error-standard-dcviation = 4%. 

With SU (4) symmetry, the picture is different. We 
have a numerical result that is reliably interpolated via 



9Dpd{s) 



5.05 - 4.26s 
1 + 0.36s - 0.060s2' 



(17) 
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FIG. 1. Upper panel - Dimensionless couplings: gopD (solid 
curve); QkpK (dashed curve); and (?,rp7r (dotted curve) - 
all computed as a function of the p-meson's off-shell four- 
momentum-squared, with the pseudoscalar mesons on-shell. 
Recall that with our Euclidean metric, > is spacelike. 
Lower panel - Ratios of couplings; gKpK / gopD (solid curve); 
and QKpK / Qirpw (dashed curve). In the case of exact SU{4) 
symmetry, these ratios take the values, respectively, 1 (dot- 
dashed line) and (1/2) (dotted line). The vertical dotted line 
marks the p-meson's on-shell point in both panels. (N.B. In 
GeV: mo = 1.865, rUp = 0.77, mn = 0.494, = 0.138.) 



A monopole parametrization at spacelike-P^, with mass- 
scale KopD = 0.69 GeV, provides a fit with relative- 
error-standard-deviation = 5%. Our computed value 
gopoiO) — 5.05 is 75% larger than an estimate ob- 
tained using QCD sum rules (3.0 ± 0.02 Q) and 100% 
larg er than a vector-meson-dominance estimate (2.52 
[36|). Moreover, if SU{4:) symmetry were exact, then 
gopD = gKpK = .97rp7r/2, but it is plain from Eq. 
that gKpKiO) = 0.92, a result which exposes a symmetry 
violation of 440% at = 0. Furthermore, on the entire 
domain € [— m^,m^], the symmetry-based expecta- 
tion gopD ~ gxpK is always violated, at a level of be- 
tween 360-440%. The second identity, gopD ~ <77rpir/2, 
is violated at the level of 320-540%. (N.B. In connec- 
tion with heavy-quark symmetry, corrections of this order 
have also been encountered c — )- c? transitions [T3|.) 

These conclusions are dramatic, so it is important to 
explain why we judge them to be robust. The computa- 



tions of gT^p-K and gxpK are considered reliable because 
we can smoothly take the limit s-quark— > w-quark and 
thereby recover a unique function that agrees with earlier 
computations by other groups. 



This leaves the possibility of uncertainties connected 
with Sc{k), Eq. HI]); rD(fc;P), Eq. dH]); and the mo- 
mentum partitioning parameter, wj. To explore sensi- 
tivity to the c-quark propagator we used an even sim- 
pler, non-confining constituent-like form; viz., Sc{k) = 
\/ {i-y-k+Mc). The effect at spacelike-P^ is modest. How- 
ever, the impact is large at timclike-P^ because there- 
upon the p-meson momentum-squared begins to explore 
a neighborhood of the spurious pole in Sc{k). Thus, 
the simpler propagator serves to increase the violation 
of SU{4:) symmetry. Regarding TD{k;P), uncertainty 
is implicit in the value of ojd ~ 1.63 ± 0.10 GeV, con- 
strained by the weak decay constant = 206± 9MeV 
[Slj . However, variations of even 20% in wd have no ma- 
terial impact on our results. Connected with that, a 20% 
change in produces only a 4% variation in via the 
fit to fD+ , hence any possibility of an effect from can 
be discounted owing to the previous consideration. 
IV. Discussion. Predictions for bound-states and res- 
onances derived from meson-exchange models are sen- 
sitive to the values of couplings in their Lagrangians. 
In these non-relativistic models the couplings are com- 
monly fixed to reproduce some known experimental data, 
e.g. the scattering length of a physical system. The 
most prominent such coupling, namely g-j^N, has long 
been used in nucleon-nucleon potentials and serves to 
define the strength of the pion's coupling to a nucleon. 
It also determines the scale of the long-range force in 
the nucleon-nucleon interaction and associated scattering 
cross sections. Analogously, the strength of the couplings 
D'-*^Dtt, i:)(*)D(*)p between D mesons and a light pion or 
p-meson plays a crucial role in the formation of charmed- 
nuclei. However, whereas g^^N can be extracted from ttA^- 
scattering data [37|, no such information is available for 
charmed-meson interactions with nucleons. 

In our approach, which is based on an internally con- 
sistent use of impulse approximation and unifies the de- 
scription of light- and heavy-mesons, we compute these 
couplings from the transition amplitude between two D 
mesons and an off-shell light meson. We find that SU (4) 
symmetry is a very poor guide to the couplings. On the 
other hand, in relation to such models it provides a con- 
structive suggestion that one might reasonably employ 



TTtME / 

^DpD\ 



ME 

yopD 



A ME 2 
^^DpD 



A ME 2 
^^DpD 



(18) 



with 



pD 



A ME 

^^DpD 



0.7 GeV, to describe DD scat- 



tering via p{q)-mesoTL exchange. 

This might be compared with the parametrization 



^Dp_d(I'?I^) — gopD 2 



aH2 
^^DpD 



DpD 



\q\ 



2 ' 



(19) 



^^DpD 



1.4 GeV, 

9 DpD ~ 2, based on the notion of 
SU{A) symmetry, which our analysis has discredited. 
The coupling in Eq. ((T^ is smaller than that in Eq. ^TE\\ 
but the evolution is harder. These effects cancel to some 
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degree, but here the magnitudes are such that our result, 
Eq. (jl8p , provides an integrated interaction 

Vo^l (20) 

that is roughly 40% greater. (N.B. If gg^^, 2.6 « 
{l/2)g^fjj, then Vq" « Vq^'^.) By the same measure, our 
DpD interaction is 20% stronger than that in Ref. [loj . 
which uses A|^^^ = 1.14GeV, gv =5.8 and hence 

glpD = 0.9gv[l - my kilo] = 2.85 . (21) 

Whilst our results argue against hard form factors, the in- 
teraction enhancement they produce is abundantly clear. 
Notably, a large value for the interaction strength entails 
an inflated cross-section in DN scattering. In particular, 
in the meson-exchange model of Ref. Q (single- meson ex- 
change version), the 1 = 1 DN cross-section is inflated 



by a factor of ^ 5, when using the our result, Eq. (|18p . 
for uj and p, instead of Eq. (jl9p . Hence, implementation 
of our results could have material consequences on, e.g., 
the possibility for formation of charmed-resonances or - 
bound-states in nuclei. 
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